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Abstract 

Collective flow of protons and negative pions has been studied within the momentum region 
of 4.2 -j- 4.5 AGeV/c (E = 3.4 -j- 3.7AGeV) for different projectile-target combinations involving 
carbon and, specifically, He-C, C-C, C-Ne, C-Cu and C-Ta. The data stem from the SKM-200-GIBS 
streamer chamber and from Propane Bubble Chamber systems utilized at JINR. The directed flow 
of protons grows dramatically in the carbon region when the counterpart nucleus grows in mass 
between He and Ta. The elliptic proton flow points out of the reaction plane and also strengthens 
as system mass increases. Within the reaction plane, the negative pions flow in the same direction 
as protons for the lighter of the investigated systems, He-C, C-C and C-Ne, and in the opposite 
direction for the heavier, C-Cu and C-Ta. The Quark-Gluon String Model reproduces observed 
changes in the flow with system mass. 
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I. INTRODUCTION 



One of the central goals of the high-energy heavy-ion collision research has been the de- 
termination of properties of nuclear matter at densities high compared to that in the ground- 
state nuclei and at temperatures high compared to energies per nucleon in the ground-state 
nuclei. Among efforts towards that goal, important discoveries made have been those of 
collective flow effects in the collisions. A collective flow is a motion characterized by space- 
momentum correlations of dynamic origin. Different aspects of flow have been investigated 
both experimentally and theoretically. The flow produces asymmetries associated with the 
reaction plane, in the particle emission patterns. Theoretically, those asymmetries can be 
linked to the fundamental properties of nuclear matter and, in particular, to the equation of 

r~i 

state (EOS) Two types of asymmetries have been identified. One has been the directed 
flow in the reaction plane, associated with the matter "bouncing-off" within the hot partic- 
ipant region of overlap between colliding nuclei. The other has been the squeeze-out of the 
hot matter moving perpendicular to the reaction plane from within the participant region. 
With energy increasing into ultrarelativistic values the squeeze-out turns into an in-plane 
elliptic flow. 

By now, the collective flow effects have been investigated over a wide range of energies, 
from tens of MeV/nucleon to 200 GeV/nucleon in the center of mass. For most part the 
experiments have relied on electronic techniques. A streamer chamber served as a detector 
in the early Berkeley experiments and later in Dubna experiments. 

Regarding investigative strategy, the reaction plane is the plane within which the centers 
of initial nuclei lie, separated in transverse direction by the impact parameter b. Spatial 
asymmetry in the initial state, associated with the reaction plane, gives rise to asymme- 
tries in the particle emission patterns. Within the method of analysis of those asymmetries, 
proposed by Danielewicz and Odyniec [2|, those asymmetries are used to estimate the di- 
rection of the reaction plane and the asymmetries themselves are assessed in relation to the 
estimated reaction-plane direction. In addition, for the analysis of asymmetries, a Fourier 
decomposition of azimuthal particle distributions has been employed relative to the 
reaction plane, in the context of particle azimuthal correlations. 

Chronologically, the collective flow of charged particles has been first observed experimen- 
tally at the Bevalac by the Plastic Ball {4], a, f| and Streamer Chamber |3] collaborations. 
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The flow continued to be explored at Berkeley and at GSI 
H, Q, El and at CERN/SPS 



AGS 
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131 ] and further at 



211 ] accelerators. The first results from 



Relativistic Heavy Ion Collider (RHIC) of BNL have been those on elliptic flow of charged 
particles at midrapidity of Au-Au collisions at the energy of ^/s = 130 AGeV 22]. 

At JINR in Dubna, emulsion techniques have been employed to investigate the directed 
flow of projectile spectator fragments in the interactions of O, Ne, Si, S, Kr, Au and Pb 
with Ag(Br) targets at beam momenta of 0.95, 4.5, 11.6 and 14.5AGeV/c [231 ] . Moreover, 
the flow of protons and n~ mesons has been studied at JINR in the central C-Ne and C-Cu 
collisions at beam momentum of 4.5AGeV/c 24j, 
SKM-200-GIBS collaboration, and in the semicentral C-C 
at beam momentum of 4.2AGeV/c, using the 2 m Propane Bubble Chamber. In those 
sets of flow investigations, on one hand the method of Danielewicz and Odyniec has been 



251 ] . using the streamer chamber of the 



26 



27j and C-Ta [28| collisions 



employed 



24 



distributions 



25 



26 



271 ] and, on the other, the method of Fourier expansion of azimuthal 



281 ] . The review [29(] has summarized some of the results obtained in the 



streamer and bubble chambers. 

In this paper, we concentrate on the collective behavior of protons and tt~ mesons in 
collisions where carbon is either a projectile or target and specifically He-C, C-C, C-Ne, 
C-Cu and C-Ta collisions at beam momentum of 4.2 H- 4.5 AGeV/c, measured either in the 
JINR streamer or bubble chambers. Details of the strategy of quantifying the collective flow 
effects differ from what the Tbilisi group had employed before 29]. 



II. EXPERIMENTAL DATA 



The data combine results obtained within the SKM-200-GIBS set-up and within the 2 m 
Propane Bubble Chamber of JINR. 

The SKM-200-GIBS setup consists of a 2 m streamer chamber placed in the magnetic 
field of 0.8 T and of a triggering system. The streamer chamber has been exposed to a beam 
of C nuclei accelerated in the JINR synchrophasotron to the momentum of 4.5AGeV/c 
(energy of 3.7 AGeV). The thickness of Cu solid target, in the form of a thin disc, was 
0.2 g/cm 2 . Neon-gas filling of the chamber also served as a nuclear target. A central trigger 
was used to select events with no charged projectile spectator fragments at p > 3GeV/c, 
within a cone of half angle of either O c h = 2.4° or 2.9°, depending on the run. The trigger 
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efficiency was 99% for events with a single charged particle in the cone. The ratio cr cen t /cinei 
(that characterizes the centrality of selected events) is 9±1 % for C-Ne and 21±3 % for C- 
Cu. Details of the acquisition techniques of 723 C-Ne and 663 C-Cu interaction data, and 
other experimental procedures such as e.g. dea ling with biases and corrections in the analysis 
have been presented in previous publications [30|, |3l| . Investigation of collective flow effects 



generally requires analyzing collisions event-by-event or exclusively. In this context, it has 
been important to put an effort into the identification of ir + mesons, the admixture of which 
amongst positive charged particles has been about 25 -=-27%. The identification was carried 
out on the statistical basis using two-dimensional (p , p L ) particle distributions [321 ] . It had 
been assumed, that tt~ and tt + mesons contribute to a given cell within the (p , p L ) plane 
with equal probability. As one detail, the difference in multiplicity of ir + and tc~ in each 
event was required to be less than 3. After the performed identification, the admixture of 
7r + mesons amongst the protons is estimated to be less than 5 4- 7%. 

The data on He-C, C-C and C-Ta interactions have been obtained using the 2 meter 
Propane Bubble Chamber of JINR. The bubble chamber was placed in a magnetic field 
of 1.5 T. Three Ta-plates (140 x 70 x 1) mm 3 in size and mounted into the fiducial volume 
of the chamber at a distance of 93 mm from each other, served as a nuclear target. The pro- 
cedures for separating out the He-C and C-C collisions in propane and the processing of 
the data including particle identification and corrections have been described in detail in 
Ref. [33(. The analysis produced 9737 He-C, 15962 C-C and 2469 unambiguously identified 
C-Ta inelastic collision events. From those, subsamples of semicentral events have been 
selected for the flow analysis, by requiring the minimal participant proton multiplicities of 
-Wpart > 3, iVp art > 4 and iVp art > 6 , respectively, for the He-C, C-C and C-Ta collisions. The 
target fragmentation products have been identified as those with momentum p < 0.3 GeV/c. 
The projectile fragmentation products have been defined as those characterized by the mo- 
mentum p > 3GeV/c and angle 9 < 4°. In consequence, 6400 He-C, 9500 C-C and 1620 
C-Ta semicentral collision events have separated out from the identified inelastic He-C, C-C 
and C-Ta collisions. The overall cross section ratios <T cen t /finei is estimated to be within the 
range of 15 4 50 % for the He-C, C-C and C-Ta collisions. 
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III. DIRECTED FLOW OF PROTONS 



In determining the directed transverse flow of protons, we have employed the method of 
Danielewicz and Odyniec . Most of the data below 4 A GeV in the literature have been, 
actually, analyzed following that method. The advantage of that method is that it can be 
employed even at small event statistics such as typical for film detectors. The method relies 
on summation over transverse momenta of selected particles in the events. 

The reaction plane is spanned by the impact parameter vector b and the beam axis. 
Within the transverse momentum method [2[|, the direction of b is estimated event-by-event 
in terms of the vector constructed from particle transverse momenta pf: 

n 

<3 = (i) 

where the sum extends over all protons in an event. The summation weight is Ui = 1 for 
Hi > y c + 8, Ui = —1 for yi < y c — 5 and uj,i = for y c — 5 < yi < y c + 5, y% is particle rapidity 
and y c is system cm. rapidity. Particles around the cm. rapidity, with weak correlations 
with the reaction plane, are not included in the reaction-plane determination. For the He-C 
and C-C interactions, we have used 5 = 0.2 and for the C-Ne, C-Cu and C-Ta interactions 
we have used 5 = 0.4. 

When referring a specific proton j to the reaction plane, we estimate the direction of the 
latter from the vector Q with contribution of proton j removed, 

to eliminate the correlation of the particle with itself, competing with the dynamic effect we 
are after. Projection of the transverse momentum of proton j onto the estimated reaction 
plane is 

The average in-plane momentum components (p x '(y)) can be obtained by averaging over 
events the momenta in different rapidity intervals. 

Due to finite number of particles used in constructing the vector Q in (CQ), the estimated 
reaction plane fluctuates in the azimuth around the direction of the true reaction plane. 
Because of those fluctuations, the average momenta (p x ') calculated with the estimated 
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reaction plane get reduced as compared to those for the true reaction plane (p x ): 

(p x '(y)) = (p x (y)) (cos $}. (4) 

Here, $ is the angle between the true and estimated reaction plane. The overall correction 
factor k = l/(cos$), which needs to be applied to (p x '(y)) in order to obtain (p x (y)), is 
subject to uncertainty, especially at low multiplicities. In this work, we evaluate (cos $) 



from the ratio 



3jj: 



{up") lupf-QA I ( {Q 2 - YIU (uPi) 



where n is proton multiplicity in an event. The (cos $) values obtained for different systems 
are listed in Table [B 

To test the quality of our reaction-plane determination, we have followed the proce- 
dure outlined in Ref. Specifically, we have divided randomly each event into two sub- 
events and we have constructed vectors Qi and Q2 for those sub-events. The distribution of 
reaction-plane directions from the sub-events peaks at zero degrees in relative azimuthal an- 
gle and it exhibits a width of a ~ 50-r55° depending on reacting system. The distribution of 
reaction planes determined from full events, about the true reaction-plane direction, should 
be about half as wide as the relative distribution of directions from the sub-events [2, [35I ] . 
The thus-estimated spreads of the estimated reaction-plane directions about the true plane 
direction, cr ~ cr/2 = 25 -j- 28° depending on a system, are comparable to values reported 



111 



361 ] and to the values of cr ~ 23 25° arrived at in Kr- and Au-induced collisions with 



Ag(Br) targets [23|. 



With regard to the rapidity, the results of our analysis of the He-C and C-C interactions 
are presented in the cm. system and those of the C-Ne, C-Cu and C-Ta interactions are 
presented in the lab. system. The values of proton (p x (y)) extracted from the collision 
events at 4.2 4.5 A GeV/c (3.4 -h 3.7 A GeV), corrected for (cos $) from ([S]), are shown in 
Fig. □ for He-C (panel a), C-C (b) and C-Ne (c) and in Fig. [2] for C-Cu (a) and C-Ta (b), 
respectively. From within this set, the He-C system is the lightest one within which the 
directed proton flow has ever been observed. The maximal values of (p x {y)) in the carbon 
region are observed to grow gradually and monotonically as the mass of the counterpart 
nucleus increases. Interestingly, the maximal values of (p x {y)) are always larger in the 
region of the lighter of two colliding nuclei. 



One commonly employed measure, quantifying the dependence of {p x )on rapidity, which 
makes no distinction between the two colliding nuclei, is the slope [5] of (p x (y)) at its 
midrapidity cross-over, 

(6) 



d(p x ) 



(p x )=0 



dy 

The slopes, represented in Figs. [1] and [2] by straight lines, have been obtained, as common, 
through third-order polynomial fits to midrapidity data, specifically from the rapidity in- 
tervals of -0.70 -T- 0.70 for He-C, -0.60 -=- 0.60 for C-C, 0.15 -=- 1.85 for C-Ne, 0.20 1.60 
for C-Cu and 0.10 4- 1.40 for C-Ta. The slopes, corrected for (cos$), are further provided 
in Table [Hand they are plotted against the geometric mean of projectile and target masses, 
in Fig. [31 

In Fig. [3] and Table (TJ it is seen that the proton F rises monotonically with the rise 
in the mass of the nucleus counterpart to C, from F p = 95 ± 8 MeV/c for He-C to F p = 
178 ±20 MeV/c for C-Ta. Within the observed mass range and errors, the dependence of'F p 
on the system mean geometric mass can be described as linear. If a power law forced through 
zero at zero mass were assumed for the dependence of F p on mass, though, a'fractional power 
would have resulted. A deviation from linearity is further anticipated as F p would need to 
saturate with mass in the hydrodynamic limit of large masses. 

The experimental results on proton directed flow have been compared to the predic- 
tions of the Quark-Gluon String Model (QGSM). Detailed description of that model can 



be found in 



37 



The QGSM is based on the Regge and string phenomenology of 
particle production in inelastic binary hadron collisions. The model oversimplifies nuclear 
effects in that nucleon mean-field effects are ignored as well as nucleon coalescence into 
clusters. Within QGSM, nuclear densities are used for selecting coordinates of original nu- 
cleons. This is followed by the formation of quark-gluon strings which fragment into hadrons. 
Those hadrons rescatter. In QGSM, the sole cause of sidewards flow is the hadron rescatter- 
ing. In all other models in the literature, where produced pressure comes from rescattering 
alone, the flow in the specific energy regime was found significantly underestimated, albeit 
predominantly in heavier systems. However, in collisions involving carbon studied here, 
QGSM turned out to describe the flow rather well. 



For simulating the model events, we have employed the COLLI Monte-Carlo generator 39] 
based on QGSM. To the generated events, a detector filter has been applied and, in the case 
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of C-Ne and C-Cu collisions, a trigger filter. In mimicking, in particular, the deterioration of 
experimental efficiency for registering vertical tracks, protons characterized by deep angles 
greater than 60° have been excluded from an analysis. In the past, it has been found that 
the applied filters selected peaked distributions of impact parameters for the collision events, 
characterized by the average 6-values of: (b) = 2.80 fm and (b) = 2.65 fm for semicentral 



He-C and C-C events 
and C-Cu 



24 



271 ]. respectively, (b) = 2.20 fm and (b) = 2.75 fm for central C-Ne 



251 ]. respectively, and (b) = 6.54 fm for C-Ta semicentral. For the present 



study, all events have been generated at fixed 6-values equal to the above averages, with 
event numbers given in Table [B Subsamples of those events have been chosen for flow 
analysis following the same multiplicity criteria as for the data. The proton flow results 
from QGSM are superimposed on data in Figs. [Hand [2J As can be seen, the model describes 
the data there rather well. The mass dependence, as seen in Fig. [3] is, in particular, fairly 
well reproduced. 



IV. PROTON ELLIPTIC FLOW 



Another flow identified in the literature, persisting even around midrapidity where the 
directed flow is suppressed, has been the elliptic flow. The elliptic flow is associated with 
the shape of participant region of overlap between the two nuclei, elliptic in the directions 
transverse to beam axis. Gradients within the elliptic participant region are stronger in 
the direction of the reaction plane than out of the plane. That anisotropy favors overall 
expansion of matter in the direction of the plane. However, at sufficiently low energies, 
when cold spectator matter lingers in the vicinity of an expanding participant matter, the 
participant expansion gets shadowed in the direction of the reaction plane and the so-called 
squeeze-out, emission out of the reaction plane dominates. Details depend on beam energy 
and on the pace of expansion of the participant matter. The energy region for our data is 
particularly interesting because it is in the vicinity of transition from squeeze-out to in-plane 
elliptic flow. Given that different particles may be emerge at different times, different elliptic 
flow directions can principally result for different particles. 

Azimuthal distributions of protons, relative to the reaction plane direction estimated 
with (j2j), are shown, respectively, in Fig. H] for He-C (a), C-C (b) and C-Ne (c) collisions 
and in Fig. [5] for C-Cu (a) and C-Ta (b) collisions. The azimuthal angle is one for which 
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cos ipj = pf/pf- The distributions exhibit a clear modulation, with maxima occurring 
around 90° and 270°, i.e. out of the reaction plane, and minima around 0° and 180°, i.e. in 
the reaction plane. This is the pattern such as observed at lower rather than higher energies. 

Notably, the fluctuations of the true reaction plane about the estimated one tend to flatten 
the azimuthal distributions of the particles relative to the plane. To quantify further the 
observed modulation, we fit the azimuthal distributions with the Fourier cosine-expansion 
(given the system invariance under reflections with respect to the reaction plane): 

diV 

— — = a (1 + a'i cos <p + a 2 cos 2<p) . (7) 
dip 

The squeeze-out signature is the negative value of the coefficient a' 2 . Compared to the 

;h a distribution relative to the true reaction plane, that coefficient 



coefficient a>2 associated wit 
is reduced, though, by 0, 4l| 



a' 2 = a 2 (cos 2$) , (8) 



in an analogy to the reduction in (j3J). The reduction coefficient for (jSJ) may be estimated 
from 

, 9(b , \((p x r - (p y, ) 2 )\ (Q , 

(cos 2$) = -rj-. — — — -rr-j- , (9) 

\((p x ) 2 - (P y ) 2 )\ 

where the numerator and denominator on the r.h.s. are, respectively, obtained from 

m 2 - ( P y T) = (2 ( ^) 2 - (pfr) , (io) 



and 



K(p*) 2 - (p y f)\ = \l {2T '' T {n Pf PiY) ■ (") 



In the above, the transverse tensor T is 



n , 

£(*>?: 

i=l ^ 



T^ = yjpt^-^fS^) , a = x,y, (12) 



and 

y 

rp . Tp _ rpa/3 rpa/3 _ ^rpxx^l _|_ (rpyy\jl _|_ ^rpxy^ 

a,/3=x 

Use of Eq. for a reaction requires that some elliptic anisotropy is present in the 
analyzed system to start with. Values of (cos 2$) estimated using Eq. flU]) for our systems 
fall within the range of 0.520 -j- 0.582, depending on the system, and are listed in Table [B 
The fits to the azimuthal distributions with Eq. ([7j) are illustrated with solid lines in Figs. 0] 
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and [51 The elliptical modulation parameters, corrected according to Eq. (JS]), from the 
fits made under different cuts to analyzed particles, are provided in Table [TTJ The elliptic 
anisotropy, quantified in terms of the a 2 -coeflicient, strengthens as the mass of a nucleus 
counterpart to carbon increases and as the transverse momentum increases. The ratio of 
out-of-plane to in-plane emission probabilities can be estimated as 

1 — 02 



R 



1 + a 2 



(14) 



with values also provided in Table [TT1 Regarding kinematic region, the elliptic anisotropy is 
seen to strengthen in Table [III as the cut around midrapidity region narrows. 

In the past, the elliptic flow has been investigated by different groups in different systems 



for different particles. Thus, the Plastic Ball 



42 



431 ] and FOPI [44J groups have investigated 



the flow in Ca+Ca, Nb+Nb, Ni-Ni, Xe-Csl, Au-Au collisions from 0.15 to 1.0 GeV/nucleon, 
for protons, light fragments and 7r =1= mesons. The KaoS 45J group has investigated the 
flow in Au+Au and Bi+Bi at 0.4-1.0 GeV/nucleon, for protons, light fragments, pions and 
kaons. Regarding theoretical descriptions, those lower-energy data have been, in particular, 



successfully described by the IQMD model 



46[ . Our results, such as of strengthening of the 



flow with increase in system mass or transverse momentum, confirm the earlier experimental 

nn 

findings, some from lower and some from higher [40J, |41[ energies. 

Our results on proton elliptic flow have been further compared to those from the QGSM 
model of the reactions. Proton distributions relative to the reaction plane estimated as in the 
experiments, normalized as in the experiments, are overlaid over the data in Figs. H] and [51 
The (^-coefficient values from the QGSM model, for distributions relative to the known 
reaction plane, are further given in Table [III As may be seen, the calculations compare 
favorably to data. On the other hand, in our testing this model appears to overestimate 
the strength of elliptic flow in the heavier Au + Au system, when the higher-energy E-895 



data 



40 



411 ] are interpolated across our energy domain. 



At higher e nerg ies, the el 



average cosine 



40 



41 



47|, 



48 



i pti c flow results have been mostly expressed in terms of the 
49] 

w 2 = (cos20), (15) 



where is the angle relative to the true reaction plane. Given the Fourier expansion (JTj), 
the average v 2 differs from the coefficient a 2 by just a factor: 

v 2 = a 2 /2 . (16) 
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Most thoroughly and over the widest energy range, the elliptic flow has been studied 
for protons in Au+Au collisions. Figure [6] compares our results to the systematics for 
Au+Au. The flow generally changes sign just above our bombarding energy It is likely, in 
the view of our results, though, that the transition energy depends on the colliding system. 
It should be mentioned that in the Au+Au measurements, due in particular to higher 
particle multiplicity, it has been possible to achieve a more narrow and more central relative 
impact parameter selection than for our systems. Determination of excitation functions 
for different-size systems could facilitate separation of the effect of transport properties of 
nuclear medium and of the equation of state and facilitate localization of the transition to 



quark-gluon plasma [501 ] . 



V. DIRECTED FLOW OF NEGATIVE PIONS 

Pions, being copious in relativistic heavy-ion collisions, can probe the reaction dynamics 
independently of nucleons. Of interest, in particular, is the relation between the directed 
flow of nucleons and pions. 

Historically, the pattern of pion emission relative to the reaction plane has been first 



studied at the _Bevalac by the Streamer Chamber group 5l|, |52|] and later by the EOS 



collaboration 



531 ] at Saturne by the Diogene group [5J]. The Diogene group investigated, 
in particular, collisions similar to ours where one nucleus, in their case the projectile Ne, 
was fixed while the counterpart target nuclei were varied. In the investigations at Bevalac 
and Saturne, projection of pion transverse momentum onto the reaction plane has been 
examined. 

For our systems, we have investigated the directed flow of negative pions. The average 
component of pion transverse momentum in the reaction plane, evaluated using Eqs. (jl]) 
and ([5]), is shown in Figs. [7] and El for He-C, C-C, C-Ne, C-Cu and C-Ta collisions, respec- 
tively. The values of flow parameter F at cross-over, from the fits to data illustrated in 
Figs. [7J and [BJ are given in Table [I] and they are further plotted against the geometric mean 
of colliding masses in Fig [31 As may observed by examining the Table and by comparing 
Figs. [TJ and M with Figs. [I] and El the situation with the directed flow of pions is different 
than with the flow of protons. The maximal transverse momentum and the flow parame- 
ter are smaller for pions than for protons. Further, while the average in-plane transverse 
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momentum of protons gradually increases in the C region as the mass of the counterpart 
nucleus increases from He to Ta, for pions the momentum first decreases in magnitude and 
then changes sign. Also, the flow parameter changes sign. 

The reduced magnitude of avera ge p ion momentum component compared to protons ha s 
been seen before at Bevalac 



54|, GSI-SIS Q and CERN-SPS 



HQ- 



, Saturne 

From the thus-far investigated systems in the literature, the maximal in-plane momentum 
and flow parameter for pions appear to be the largest in magnitude relative to protons for 
our C-Ta system, where F 71 = —74 ± 7MeV/c. The direction of pion flow opposite to 



oroton flow, termed antiflow, has been seen before in either asymmetric [54j or symmetric 



20 



531 ] systems. However, we are unaware of an observation of the pion antiflow in 



strongly asymmetric systems, such as our C-Cu and C-Ta, where the pion in-plane mo- 
mentum would be simultaneously changing sign as a function of rapidity. The cross-over in 
in-plane momentum for pions is pushed towards C or nucleon-nucleon midrapidity compared 
to protons in our strongly asymmetric C-Cu and C-Ta systems. 

If only the collective motion played a role in giving rise to the finite average in-plane 
particle momenta, the pattern of (p x ) vs rapidity should be similar for protons and pions, 
with (p x ) /m being of comparable magnitude for the particles 55J. This is actually what is 
observed for our lighter systems, He-C, C-C and C-Ne. The comparable values of particle 
(p x )/m imply, though, small values of pion (jf) due to collective effects, that can be then 



easily modified by an shadowing effects [55 



56, 



57| . The latter effects can become pro- 



nounced in more peripheral collisions of heavier systems. Successful restricting of relative 
centrality in a collision relying on multiplicity gets difficult in highly asymmetric systems. 
With this, the conditions become ripe for shadowing in our C-Cu and C-Ta systems. Quanti- 
tatively, the shadowing can produce in-plane transverse momentum components comparable 
to the momenta itself and, thus, much larger than components due to collective motion for 



pions 



]. As can be seen from the QGSM results overlaid over the data in Figs. [7] and [HI 
the model reproduces the pion flow-antiflow effects found in the experiment. 

Given the past and present results for asymmetric systems, it could be useful to study 
highly asymmetric systems at higher statistics allowing for finer separations of centrality, 
to illuminate when antiflow in the region of the lighter nucleus might change into flow, at 
antiflow persisting in the rapidity region of heavier nucleus. 
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VI. CONCLUSIONS 



Characteristics of transverse collective flow for protons and n~ mesons emitted from He- 
C, C-C, C-Ne, C-Cu and C-Ta reactions at momentum 4.2 -e- 4.5AGeV/c (energy 3.4 -e- 
3.7AGeV) have been determined. The He-C system is the lightest in which either the 
directed or elliptic flow has been ever detected. 

In determining the directed flow of protons and pions, the transverse momentum tech- 
nique of Danielewicz and Odyniec has been employed. Proton in-plane momentum appears 
always to maximize in the region of the lighter of the nuclei. In the C rapidity-region, 
maximal value of proton average in-plane momentum increases dramatically as mass of the 
counterpart nucleus increases. Also, the flow parameter at rapidity cross-over increases, 
from F p = 95 ± 8 MeV/c for He-C to F p = 178 ± 20 MeV/c for C-Ta. 

Negative pions exhibit directed flow consistent with that for protons in He-C, C-C and 
C-Ne collisions. In those collisions, the relation between the momenta is consistent with 
that expected for pure effects of collective motion, i.e. {p x ) w > {m % /m p ) (p x ) p and F w > 
(m^/rrip) F' p . On the other hand, for the heavier systems, C-Cu and C-Ta, pion flow turns 
into antiflow, with pion average in-plane momenta becoming opposite to those for protons 
both in the C and counterpart nucleus rapidity regions. The maximal pion momenta achieve 
larger values than in the lighter systems, likely due to shadowing, with magnitude of pion 
flow parameter reaching 40% of that for protons, at F n = —74 ± 7MeV/c for C-Ta. 

In assessing the proton elliptic flow, we have combined Fourier decomposition of proton 
azimuthal distributions with a transverse tensor method. In all investigated systems at our 
energies, the elliptic flow points out of the reaction plane, as is characteristic for reactions 
at lower beam energies. The strength of the elliptic flow increases as mass of the system 
increases, from a 2 = -0.039 ± 0.008 for He-C to a 2 = -0.093 ± 0.009 for C-Ta. 

The flow measurements have been compared to the QGSM model of collisions. The model 
adequately describes flow in the measured systems. 
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TABLE I: Characteristics of the measured collision events and those simulated within the 
Quark-Gluon String Model (QGSM), including event number N prior to multiplicity cut. 





He-C 


C-C 


C-Ne 


C-Cu 


C-Ta 


-^cxpt 


9737 


15962 


723 


663 


2469 


^QGSM 


25000 


50000 


16337 


5137 


10000 


(cos <J>) 


0.866 


0.893 


0.848 


0.784 


0.720 


(cos 2$) 


0.563 


0.582 


0.531 


0.559 


0.520 


Fl pt (MeV/c) 
(^-corrected) 


95 ± 8 


115 ± 11 


123 ± 12 


143 ± 15 


178 ± 20 


^qgsm (MeV/c) 


93 ± 6 


111 ± 8 


125 ± 8 


138 ± 11 


183 ± 12 


F- pt (MeV/c) 
(fc-corrected) 


17 ± 3 


19 ± 3 


24 ± 5 


-43 ± 6 


-74 ± 7 


^qgsm (MeV/c) 


18 ± 3 


18 ± 3 


23 ± 3 


-43 ± 4 


-72 ± 4 
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TABLE II: Characteristics of proton elliptic flow. 



System 


Applied Cut 


expt 
^2 


R 

■*■ ^expt 


QGSM 




— 1 3 <C <C 1 3 


-0.039±0.008 


1.081±0.025 


-0.038±0.006 


He-C 


— 1.3 < lir-rr, < 1.3; 

p x > 0.2 GeV/c 


-0.049±0.010 


1.103±0.040 


-0.047±0.007 




-0.4 < < 2.1 

^ > 0.3 GeV/c 


-0.062±0.011 


1.132±0.024 


-0.056±0.007 




-1 3 < u„™ < 1 3 

_ yam — 


-0.048±0.005 


1.101±0.025 


-0.048±0.004 


C-C 


— 1.3 < Vnrr> < 1.3: 

p x > 0.2 GeV/c 


-0.064±0.006 


1.137±0.040 


-0.064±0.004 




-0.4 < Vrrr, < 2.1 

p 1 - > 0.3 GeV/c 


-0.076±0.007 


1.164±0.032 


-0.075±0.004 




-0 4 < V/„f, < 2 1 


-0.052±0.016 


1.110±0.040 


-0.055±0.004 


C-Ne 


-0.4 < y/„h < 2.1; 
p 1 - > 0.2 GeV/c 


-0.067±0.020 


1.144±0.039 


-0.069±0.004 




-0.4 < Vlnh < 2.1 

p 1 - > 0.3 GeV/c 


-0.080±0.020 


1.174±0.049 


-0.081±0.004 




-0 4 < < 2 1 

v ' • ^ _ yiao — 


-0.072±0.013 


1.155±0.038 


-0.070±0.005 


C-Cu 


-0.4 < vinh < 2.1; 
p 1 - > 0.2 GeV/c 


-0.077±0.015 


1.167±0.047 


-0.080±0.005 




-0.4 < Vlnh < 2.1 

p x > 0.3 GeV/c 


-0.093±0.016 


1.206±0.060 


-0.094±0.005 




-0.3 < y;„ fc < 2.0 


-0.093±0.009 


1.205±0.062 


-0.095±0.004 


C-Ta 


-0.3 < y ia6 < 2.0; 
p L > 0.2 GeV/c 


-0.123±0.010 


1.280±0.065 


-0.136±0.004 




-0.3 < y lab < 2.0 
p 1 - > 0.3 GeV/c 


-0.152+0.012 


1.358±0.075 


-0.162±0.004 
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FIG. 1: Average component of proton transverse momentum in the reaction plane, as a function 
of rapidity in the cm. system of He-C (a) and C-C (b) collisions and in the lab. system of C-Ne (c) 
collisions. The data, corrected for (cos <J>) from Table HI are represented by circles. Results from 
the QGSM model are represented by crosses. Straight-line stretches represent the slope of data 
at midrapidity cross-over, obtained by fitting the data with a third-order polynomial within the 
rapidity region of -0.70 < y < 0.70 for He-C, -0.60 < y < 0.60 for C-C and 0.15 < y < 1.85 
for C-Ne. The curved lines guide the eye over data. 
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FIG. 2: Average component of proton transverse momentum in the reaction plane, as a function 
of rapidity in the lab. system of C-Cu (a) and C-Ta (b) collisions. The data, corrected for (cos <3?) 
from Table HI are represented by circles. Results from the QGSM model are represented by crosses. 
Straight-line stretches represent the slope of data at midrapidity cross-over, obtained by fitting the 
data with a third-order polynomial within the rapidity region of 0.20 < y < 1.60 for C-Cu and 
0.10 < y < 1.40 for C-Ta. The curved lines guide the eye over data. 
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FIG. 3: Dependence of the slope F Q on geometric mean of projectile and target masses (Ap * 
At) 1 ' 2 , for the systems indicated in the figure. The data are represented by circles and triangles for 
protons and pions, respectively. The lines represent linear fits to the data. Results of the QGSM 
model are represented by stars in the top and bottom part of the figure, for protons and pions, 
respectively. 
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FIG. 4: Azimuthal distribution of protons relative to the estimated reaction-plane within the cm. 
rapidity region of —1.3 < y < 1.3 in He-C (a) and C-C (b) collisions and within the lab. rapidity 
region of —0.4 < y < 2.1 in C-Ne collisions (c). The circles represent data for 20 bins of azimuthal 
angle. The lines represent fit to the data with the function dN/dtp = ao(l + a'i cos (p + a' 2 cos 2ip). 
The stars represent distribution calculated within the QGSM model, with a normalization set to 
match that of the data. 
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FIG. 5: Azimuthal distribution of protons relative to the estimated reaction-plane within the lab. 
rapidity region of —0.4 < y < 2.1 in C-Cu (a) and —0.3 < y < 2.0 in C-Ta collisions (b). The circles 
represent data for 20 bins of azimuthal angle. The lines represent fit to the data with the function 
dN/d(f = ao(l + a'! cos ip+a' 2 cos 2ip). The stars represent distribution calculated within the QGSM 
model, with normalization set to match that of the data. 
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FIG. 6: Excitation function of elliptic flow in terms of the ^-coefficient vs laboratory energy per 
nucleon. The triangles represent our results slightly offset in energy to make it easier to discern the 
results for different systems. The symbols representing results of other collaborations, for Au-Au, 
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FIG. 7: Average component of 7r _ -transverse momentum in the reaction plane, as a function of 
rapidity in the cm. system of He-C (a) and C-C (b) collisions and in the lab. system of C-Ne (c) 
collisions. The data, corrected for (cos from Table HI are represented by circles. Results from 
the QGSM model are represented by crosses. Straight-line stretches represent the slope of data 
at midrapidity cross-over, obtained by fitting the data with a third-order polynomial within the 
rapidity region of -0.90 < y < 0.90 for He-C, -0.75 < y < 0.75 for C-C and 0.15 < y < 1.20 
for C-Ne. The curved lines guide the eye over data. 
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FIG. 8: Average component of 7r~-transverse momentum in the reaction plane, as a function of 
rapidity in the lab. system of C-Cu (a) and C-Ta (b) collisions. The data, corrected for (cos <3?) 
from Table HI are represented by circles. Results from the QGSM model are represented by crosses. 
Straight-line stretches represent the slope of data at midrapidity cross-over, obtained by fitting the 
data with a third-order polynomial within the rapidity region of 0.15 < y < 1.60 for C-Cu and 
0.0 < y < 1.25 for C-Ta. The curved lines guide the eye over data. 
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